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ABSTRACT 

We analyze optical and near-infrared data of a sample of 11 barred spiral galaxies, in order 
to establish a connection between star formation and bar/spiral dynamics. We find that 22 
regions located in the bars, and 20 regions in the spiral arms beyond the end of the bar present 
azimuthal color/age gradients that may be attributed to star formation triggering. Assuming a 
circular motion dynamic model, we compare the observed age gradient candidates with stellar 
populations synthesis models. A link can then be established with the disk dynamics that allows 
us to obtain parameters like the pattern speed of the bar or spiral, as well as the positions of 
resonance radii. We subsequently compare the derived pattern speeds with those expected from 
theoretical and observational results in the literature (e.g., bars ending near corotation). We 
find a tendency to overestimate bar pattern speeds derived from color gradients in the bar at 
small radii, away from corotation; this trend can be attributed to non-circular motions of the 
young stars born in the bar region. In spiral regions, we find that ~ 50% of the color gradient 
candidates are "inverse", i.e., with the direction of stellar aging contrary to that of rotation. The 
other half of the gradients found in spiral arms have stellar ages that increase in the same sense 
as rotation. Of the 9 objects with gradients in both bars and spirals, 6 (67%) appear to have a 
bar and a spiral with similar $7^, while 3 (33%) do not. 

Subject headings: galaxies: kinematics and dynamics — galaxies: star formation — galaxies: spiral — 
galaxies: photometry — galaxies: stellar content — gala:xies: structure 

1. Introduction. 



^Visiting astronomer at Kitt Peak National Observa- 
tory, National Optical Astronomy Observatory, which is Since t he earliest A''-body simulations (e.g.. 



operated by the Association of Universities for Research |Hohllll97ll ) established that bars can form in Spi- 



in Astronomy (AURA), under cooperative agreement with ^j^j^ ^^^^ theoretical work has been un 

the National science foundation. 
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Observatory, National Optical Astronomy Observatory evolution. Theoretical models show that bars Can 

which is operated by the AURA, under contract with the 
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be formed a s part of the natura l development of 
the system ( Athanassoulal lioooh . On the clas- 
sical view, bars are mainly supported by elon- 
gated orbits_(caJledxi_orbit along the bar major 
axis ( Skokos et al. 20021 and references therein). 
Chaotic orbits, however, ca n, also support bars 



in disk galaxies (se e, e.g., Voglis et aD l2006at 



Harsoula et al.ll2011l ). Given that the stars that 



make up bars remain most of the time within 
them, bars, unlike spiral arms, are not density 
waves (Sparke & Gallagher 2007; see also 15.31) . 
Bars, however, can also be considered as long-lived 
waves and "normal" m odes in the disk, possi bly 
driving spiral structure ( Buta fc Combes! 19961 ). 

It has been predicted that single large-scale 
( "fast" ) bars commonly end inside corotation (CR 



Contopoulos '1980*;' Co ntopoulos fc P apavan nopoulosI 
1980 : ScUwoo d .1981.), also called Lagran ge ra- 
dius (see, e.g., Sellwood fc Wilkinson 19931) . This 
prediction has al so been co rroborated by observ a- 
tions (e.g...Merrifield fc Kuiiken 1995: Gerssen et~al] 
19991: iDebattista et all 120021: lAguerri et ahl 120031: 



Fathi et al.l 120091 : ICabbasov et all 120091 : ICor'siiiil 
201flf ). From these studies, nonetheless, it can also 
be inferred that the bar's corotation radiuf0 lies 
between 1 and 1.4 times the bar's semimaior axis. 



and not exactly where the bar ends ( Aguerri et al 
20031 : iButa fc Zhanj2009tlElmegreeiJll996[ ). This 



may explain why some galaxies (e.g., NGC 1300, 
NGC 1365, NGC 5236) show dust lanes on the 
inside of the spiral arms and HII regions on the 
outside. // bar and arms have the same pattern 
speed, this makes sense only whi le spirals lie within 
corotation ( Roberts et al. 1979t ). A crossing of the 
dust lanes from the inside to the outside of the 
arms, marking the corotation position, or an arm 
bifurcatiord are also sometimes observed. 

Another theory postulates that "slow bars"[l 
actua lly end near the inner Lind blad resonance 
(ILR: IPasha fc PolvachenkolfiggJ )- and that they 
are formed by the alignment of elon gated and os- 
cillat ing orbits along the poten ti al (iLynden- Bdl 
19791) . ICombes fc ElmegreenI (11993^ 



other hand. 



on the 

based on the results of numer- 



^ ICeverino fc KlvpinI l l2007t) show that corotation particles 

are actually located in a "wide" ring. 
^Bifurcations are not a unique signature of CR, and can also 

be triggered by other resonances (e.g., 4:1). 
•^Those for which the CR ra dius is larger than 1 .4 times the 

bar's semimajor axis, (e.g.. lAeuerri et al.ll2003) . 



ical models (see also iRautiainen et al. 2005f) . 
propose that bars in early Hubble type galax- 
ies (those with high bulge to disk mass ratio) 
end near CR, whereas those in late type galax- 
ies (with low bulge to disk mass ratio) may 
end near the ILR, where the sp iral arms be- 
gin ()Elmegreen fc Elmegreenlll989l) . Nevertheless, 
"slow" rotating bars are not favored by the com- 



paris o n of s i niulations and observations (ISchwarz 
'1984 Il985t IWeiner et all l200lL iDubinski et al 



120091 ) ■ probably implying that disks in barred 
galax ie s are maximal ([Debattista fc Sellwood 
20001 ). iButa fc Zha"^ (120091 ) apply the potential 



density phase-shift method (jZhang fc Butal 120071 ) 
to a near-infrared (NIR) subsample of the Ohio 
State University Bright G alaxy Survey (OS- 
UBGS, lEskridge et aLll2002l ). and also fail to find 
evidence of the existence of "slow" bars. 



Sellwood fc Sparkd (|1988l ) propose that mul- 
tiple pattern speeds may b e common (see also 



Sellwood fc WilkinsoJ Il993l: IPatsis et al.l |2009[) . 



such that for some barred galaxies two corotation 
radii, one belonging to the bar and the other to 
the spiral arms, may occur simultaneously in the 
disk. This would imply that spiral arms are con- 
nected to bars only transitorily. The possibility 
also exists that the bar and bulge mask the spiral 
arms at small radii; the arms would only start to 
be seen near the bar's end, creating the impression 
of a physical connection where there is none. Un- 
sharp masking and/or Fourier-based methods can 
be used to search for spiral and bar perturbations 
in the inner parts of galaxies. However, ring-like 
features at the end of the bar may interfere with 
the unambiguous identification of a physical con- 
nection between bar and arms in some objects. 



Rautiainen fc Salo (I1999I ) perfor m simulations 



that c onfirm the scenario advanced bv lSellwood fc Sparke 
(jl988l ). They also find cases, however, where 
the bar and spiral rotate with the same pattern 
speed, as well as other cases where the bar and 
spiral patte rns have different speeds but coupled 



reson ances (jTagger et al.l Il987t iMasset fc Tagger 



19971 ). In these cases, the bar's CR may overlap 



with the spiral 4:1 resonance or ILR, and the bar 
and spiral arms are indeed connected. Various 
simultaneo us spiral modes are a lso possible. Ac- 
cording to Buta fc Zhangj (|2009t ). barred galaxies 
can have "fast" bars (CR at 1 to 1.4 times the 
bar extent), and spirals rotating at the same or 
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with a different pattern speed. They also consider 
"super-fast" bars, for which the bar's CR is well 
inside the bar's end and spirals are decoupled, al- 
though so far no theoretical models (e.g., n-body 
or response) have produced a consistent system 
with CR within the bar. 

Observations also indicate that around 20% of 
early spi ral galaxies are double-bar or nested bar 
systems l|Laine et al.ll2002t lErwinll2Q04 12009[ ). In 
these systems, an outer or primary bar harbors an 
inner or secondary bar in the nucleus. The inner 
bar rotation is independent from that of the outer 
bar. There have also been suggestions of triple-b ar 
systems (e.g. NGC 2681. lErwin fc Snarkel[l999h . 



star formation activity in the bar and cen tral re- 
gions (see alsolGarcfa-Barret o et al. 19961 ). 



1.1. Star formation 

It is commonly observed that giant HII region 
complexes occur at the end o f the bar, where the 
spiral arms seem to or iginate ([Roberts et al. I ll979t 
Kennev fc LordI Il991r ). However, not all barred 
spirals show this feature clearly, and they present 
Ha emission in other regions as well, e . g., along 
the bar (see Garcfa-Barreto et al.l 1996[ ) . PhillipsI 
(|l993l 119961 ) proposes that barred galaxies can be 
classified according to their star formation proper- 
ties. Galaxies with Hubble type SBb and "flat"El 
bars have scant star formation in the bar region 
and a concentration of HII regions in inner ringsH 
Star formation activity can be very high in the cir- 
cumnuclear region. The second group corresponds 
to Hubble type SBc and "exponential" bar type. 
This group has luminous HII regions in the bars, 
poorly defined rings, and less star formation in the 
circumnuclear region. In both SBb and SBc galax- 
ies, the star formation rates per unit area seem to 
be enhanced in the region where the bar joins the 
spiral a r ms, m ore noticeably in the SBb group. 
PhillipsI (jl996l ) points out that SBa galaxies tend 
to show star formation in ring structures, but no 



■^The light profile of flat bars decreases with radius more 
slowly than the spiral arm profile. Conversely, the radial 
scale length of exponent ial bars is the same or shorter t han 
that of the spiral arms llEImeereen fc ElmegreenlllQSSl ) 

^ Inner rings are commonly found at the end of the bar (e.g., 
lAthanassoula et al.ll2009a^ . Barred spirals may also show 
rings in the nucleus and near the outer Lindblad res- 
onance dSellwood fc Wilkinsoiil Il993t) . Rings are com- 
monly associated both with resonances llSchward llQSll. 
11984 iBvrd et al. II199 4Y and with active and recent star 
formation l lButa fc Combesi,1996i) . 



According to .Athanassoulai (|1992l ) , gas density 
enhancements perpendicular and at the end of the 
bar are often seen in simulations. She also no- 
tices that, for bars with "straight" dust lanes, star 
formation is difficult along the bar because of con- 
siderable shear (the situation may be different for 
curved lanes/shocks). These "stra ight" dust lanes 



occur in strongly barred g alaxies (IConieron et al 
2009t ). Regardless of this, ISheth et al.l (|200d l) 



ar- 
gue that some stars may be formed between the 
bar's end and the circumnuclear region. In their 
model, stars are born in "dust spurs" on the trail- 
ing side of the b ar th at feed the main dust lanes. 
Zurita fc Perez ( 20081 ). based on iJ^, optical, and 
NIR observations, find evidence to support this 
model in the barred spiral NGC 1530. They sug- 
gest that stars are formed in the dust spurs and 
migrate across the b ar to its le ading side (see also 
lAsif et al. 2005: E lmegreen et a l. 2009). 

The existence of phase shifts between the 
emission, the old stellar component of the bar, 
an d the CO emissio n alon g it has been notice d 
by iMartin fc Friedlil (Il997l) . IVerlev et al.l (|2007l ). 
and lSheth et al.l (|2002l ). The leads both the 
CO and old stars, although no systematic pattern 
with galactocentric radius is measured. These ob- 
servations support the idea that star formation can 
be triggered in the bar region. 

If spiral arms in barred galaxies trigger star 
formation, then azimuthal color gradients must 



be observed across them (e.g. Roberts _ 1969t 



Gonzalez fc Graham! 1996; Martfnez-Garcfa et al 
2009aM b[) ■ due to the aging of stars and their veloc 



ities relative to the spiral arms. In the cases where 
spirals corotate with the bar, inverse color gradi- 
ents (i.e., the direction of stellar aging is contrary 
to the sense of rotation) are expected. 

Age gradients across bars have als o been pre- 

diete d by recent n umerical simulat i ons ([Dobbs fc Pringle 
2010l ). Moreover. iMazzuca et all (|2008l) have de- 



tected azimuthal age gradients in ~ half of a sam- 
ple of 22 nuclear rings, that they speculate may 
be related to enhanced star formation in the con- 
tact points of bar and ring. In the bar region, 
the mechanisms triggering star formation may be 
different from the one that sets off star formation 
in spiral arms, due to the diverse dynamics and 
shock conditions. 
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Here, we aim to investigate azimuthal color 
(age) gradients across bars and spirals in SB (or 
strong SAB) galaxies, and the connection between 
these color gradients and bar/spiral dynamics. In 
figure m we show expected sites of color gradients; 
in this example, the bar ends near corotation (CR) 
and spiral arms corotate with the bar, which is not 
necessarily always the case. 

2. Observations 

Our sample consists of 11 barred spirals, 9 of 
them classified as SB, and two as SAB in the Third 
Reference Catalogue of B right Galaxies (RC3, 



de Vaucouleurs et al.l 119911) . Statistics show that 



^30 % of spirals are strongly barred in the opti- 
cal (de Vaucouleur^ ll963HSellwood fc Wilkinso3 



19931 ). In the NIR. lEskridge et al.l (|2000[ ) find that 
^--^56% are strongly barred (SB) and ~16%, weakly 
barred (SAB). Their study also shows that late 
type spirals (Sc-Sm) hold the same fraction of bars 
as early types (Sa-Sb). 

Optical and NIR data were acquired during 
1992-1995 with four different telescopes in three 
observatories: the Cerro Tololo Interamerican Ob- 
servatory (CTIO) 0.9-m and 1.5-m, the Lick Ob- 
servatory 1-m, and the Kitt Peak National Ob- 
servatory (KPNO) 1.3-m telescopes. Deep photo- 
metric images were taken in the optical filters g, 
r, and i, and in the NIR J, Kg or K'. Eff'ective 
wavelengths and widths of all the filters are listed 
in table [TJ the observation log for the 11 galaxies 
is shown in table [2] 

The CTIO 0.9-m optical telescope used a Tek 
10242 and a Tek 2048^ CCDs, both with a 0'.'4 
pixel^^ plate scale. The CTIO infrared obser- 
vations were performed at the 1.5-m telescope, 
with the Ohio State Infrared Imager/Spectrometer 
(OSIRIS) in February 1994, and with the CIRIM 
instrument in September 1994 and September 
1995. Both CIRIM and OSIRIS used 256^ NIC- 
M0S3 arrays; the CIRIM focus was adjusted to 
give a 1'.'16 pixel"^ plate scale, whereas OSIRIS 
was used in the /3 image mode, that provided a 
plate scale of I'.'l pixel" ^ The CCD at the Lick 
1-m telescope was a Ford 2048^, with a plate scale 
of 0'.'185 pixel"'^. For the infrared observations at 
Lick, the same telescope was fitted with the LIRC- 
2 camera; it had a 256^ NICMOS II detector, with 
a 1'.'145 pixel" 1 plate scale. The KPNO infrared 



observations were made with the IRIM camera, 
that employed a 256^ NICM0S3 array, with a 2" 
pixel" "'^ plate scale. 

2.1. Data reduction 

The data were properly r educed with the im age 
processing package IRAPO (|Todvlll986l Il993[ ) and 
fortran 77 routines. In the optical we applied over- 
scan and trimming corrections. We produced a 
combined bias frame with the "minmax" rejection 
algorithm and subtracted it from the data. Dark 
current frames were inspected, but their counts 
were found to be negligible. Twilight flats were 
compared with dome flats, and the former were 
found to be better. A "master" flat field image 
was obtained for each filter by scaling flat fields 
by their median, averaging them with a sigma 
clipping algorithm, and dividing the result by its 
mean. Objects were then divided by the appropri- 
ate "master" flat. Bad pixels masks were created 
with the assistance of both dome and twilight flat 
images. Sky subtraction was achieved by masking 
bright objects and obtaining the median of the re- 
maining pixels. For some objects the sky emission 
was not uniform and we fitted a plane, instead of a 
constant median value. To produce mosaics, indi- 
vidual images were superpixelated by a factor of 2 
in each dimension. Each superpixelated frame was 
inspected for artifacts (e.g., asteroid traces) and 
the respective pixels were masked. The images 
were then registered to the nearest pixel (i.e., half 
an original pixel), and a median mosaic was ob- 
tained. Cosmic rays were located and masked by 
comparing the median mosaic with each individual 
frame. A final mosaic was then obtained by adding 
the superpixelated, registered, clean frames. 

For the NIR data, we first corrected for the non- 
linearity of the detector. A polynomial function 
was adjusted to each pixel of dome flats of increas- 
ing exposure times. Source variability and read- 
out delay time images were used when available, 
or obtained via iterations. The flat-field correction 
was done analogously to the optical case. Sky and 
object frames were taken in an alternating fashion, 
and with no more than a few minutes difference; 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the National Science Foundation. 
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Fig. 1. — Schematic of stellar age gradients in a barred galaxy. Gradients are shown by arrows that go from 
blue (young stars) to red (old stars), across a strong bar with straight dust lanes, and through the spiral 
arms; they are produced by stars born in a shock that then drift away as they age. In this example, the bar 
ends near CR and the arms corotate with the bar; because the observed arms are completely beyond the CR 
radius, the spiral pattern overtakes the gas in the disk and the gradients in the arms are "inverse" , that is, 
the stellar aging vector is opposite to the direction of disk rotation. 
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sky exposures were centered at different positions 
with respect to the objects, and separated from 
them by ~ 10', or about twice the hnear size of 
the field-of-view. Bright objects in the sky frames 
were masked, and the four sky exposures closest 
in time to each object were median scaled, aver- 
aged with a rejection of deviant pixels, scaled to 
the mean sky value of the object, and subtracted 
from it; this process also takes care of dark cur- 
rent removal. Flat field corrections were applied 
with master flats derived from dome flats. Mosaic 
registering and cosmic ray masking were achieved 
with the same procedure used for the optical data. 



2.2. Calibration 



The optical calibration was done in the Thuan- 



Gunn system (jThuan fc GunnI 119761 : IWade et al 
19791) . The zero point of this photometric 
system is chosen such that the standard star 
BD-|-17°4708 has g = r = i = 9.5. Syn- 
thetic magnitudes were obtained for other spec- 
trophotometric stand ard^ using the spec t ral en- 
ergy d i stributions in lOke fc GunnI (1983 ). Stone 
' 1977|) . iMassev et all d 19881) . i Massev fc Gronwall 



1990 ). iHamuv et al.l (|1992D , and iHamuv et al 



19941 ) , and system response curves for each detec- 



tor /filter /te lescope/observatory combinat ion (for 
details, see iMartfnez-Garcfa et al. 2009a ). The 
best photometric galaxy frame was selected for 
each object, and the final mosaic was scaled to 
it. The NIR J, H and Kg data were calibrated 
with frames from the Two Micron All Sky Sur- 
vey ('2MASS. [Sbutskie et al.l [ l997i l2006i). The K 
data were calibrated with the ICarter fc Meadows! 
(11995 *) photometri c standard st ars, transformed to 
the CTIO system (|Carteij|l990l ). For the K' data, 



we ad opt K' = K+0.2{H-K) (|Wainscoat fc Cowid 
19921) ■ and use the pho tometric standard stars 
of Hawarden et al.l ( 200ll ). The optical data were 
degraded to the resolution of the NIR data. 



^ Feige 15, 25, 34, 56, 92, 98; Kopflf 27; LTT 377, 7987, 9239; 
EG 21; BD-(-40°4032; and Hiltner 600. 
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Table 1 

Filter char acteristics 



Filter 




FWHM 


g 


5000A 


830A 


r 


6800A 


I330A 


i 


7800A 


1420A 


J 


1.25^1111 


0.29/jm 


H 


1.65^m 


0.29/2m 


Ks 


2. 16^m 


0.33/^m 


K' 


2.11^m 


0.35/im 


K 


2.2jum 


0.41/im 
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Table 2 
OesEmATiON Log 



Object 



Filter 


Exposure(s) 


Telescope 


Date (month/year) 


Mean seeing 


g 


8100. 


Lick 1 m 


10/92, 11/92, 9/93 


2.1 


r 


7200. 




10/92, 9/93, 11/93 




i 


4200. 




10/92, 9/93 




J 


672. 


Kitt Peak 1.3 m 


11/94 




Ks 


336. 






g 


3865. 


Lick 1 m 


9/93 


1.7 


r 

i 


4695. 
3600. 




9/93, 11/93 
11/93 




J 


1380. 


Kitt Peak 1.3 m 


9/93, 11/94 




Ks 


360. 




9/93 




g 


5400. 


Lick 1 m 


3/93, 4/93, 4/94 


1.9" 


r 


7599. 




4/93, 2/94, 4/94 




i 


2100. 




2/94, 4/94 




J 


1252. 


Lick 1 m 


2/95 




J 


840. 


Kitt Peak 1.3 m 


3/94 




Ks 


540. 






g 


14646. 


Lick 1 m 


11/92, 9/93, 11/93, 10/94, 11/94 


1.6 


r 


10165. 




11/92, 9/93, 10/93, 11/93, 10/94, 11/94 




i 


6600. 




11/92, 9/93, 11/93, 10/94, 11/94 




J 


300. 


Kitt Peak 1.3 m 


9/93 




k' 


352. 


Lick 1 m 


12/94 




g 


3600. 


CTIO 0.9 m 


9/94 


1.5 


r 


5100. 








i 


3900. 








J 


810. 


CTIO 1.5 m 


9/94, 9/95 




Ks 


420. 








g 


3600. 


CTIO 0.9 m 


9/94 


1.6 


r 


3600. 








i 


3600. 








J 


300. 


CTIO 1.5 m 


9/95 




Ks 


330. 








g 


2580. 


Lick 1 m 


4/94, 11/94 


1.4" 


r 


2707. 








i 


780. 








J 


626. 


Lick 1 m 


2/95 




J 


840. 


Kitt Peak 1.3 m 


3/94 




Ks 


560. 






g 


3600. 


CTIO 0.9 m 


3/94, 3/95 


1.5" 


r 


4500. 








i 


4200. 








J 


971. 


CTIO 1.5 m 


2/94 




K 


868. 








g 


4800. 


CTIO 0.9 m 


3/94, 3/95 


1.6" 


r 


4200. 








i 


5400. 








J 


313. 


CTIO 1.5 m 


2/94 




H 


313. 








g 


10620. 


Lick 1 m 


8/92, 9/93, 10/94 


1.8" 


r 


6960. 








i 


8649. 








J 


1080. 


Kitt Peak 1.3 m 


9/93, 11/94 




Ks 


396. 








g 


4200. 


CTIO 0.9 m 


3/94, 3/95 


1.5" 


r 


4500. 








i 


4200. 








J 


2719. 


CTIO 1.5 m 


2/94 




H 


2718. 






K 


1042. 









NGC 718 , 



NGC 864 



NGC 4314 



NGC 266 , 



NGC 986 , 



NGC 7496 



NGC 5383 



NGC 4593 



NGC 3059 



NGC 7479 



NGC 3513 
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3. Data analysis 

The objects were deprojected, by first ro- 
tating the frames to align the galaxy's major 
axis with the "y" direction, and then stretch- 
ing the rotated image in the "x" direction by 
a factor 1/cosa, where a is the inclination an- 
gle. The position and inclination angles were 
taken from the literatur e, mostly from Hyper- 
leda (|Paturel et al. 2003 ) and the RC3 (see ta- 
ble |3l). The bar and spiral arm perturbations 
were traced in the NIR bands (mainly Kg, K', 
and K, but other when specified), assuming that 
young stars do not contribute t o the observed ra- 
diation fe.g-. lRix fc Riekelll993l) . However, young 
stars and clusters may contribute l ocally up to 
20% — 30% to the observed radiation (IRix fc Riekd 
19931: iGonzalez fc Grah"a3 Il996t iRhoadd Il99 



Patsis et al.l2OOlHGrosb0l et alfeOOeHCrosbol fc Dottori 
20081) . 



the Q{rJgi) "unwrapped" frame, with the aid of 
a dust lane tracer (e.g., the g — J color). The 
arms are then "straightened" , by adding a dif- 
ferent phase shift at each radius, until the arms 
appear "horizontal" . In the straightened arms, 
selected regions can be collapsed in In R to yield 
1-D plots of intensity vs. distance that can be 
compared with stellar population synthesis (SPS) 
models. 

Stellar models give Q{rJgi) as a function 
of stellar age (tagc), while observations provide 
Q{rJgi) as a function of distance. Distance d = 
is fixed at the location with the highest extinction 
(i.e., the highest g — J value) in the dust lane. 
The pattern speeds of the bar, fi^'^'', and of the 
spiral arms, ripP"'^\ are derived from candidates 
in, respectively, the bar and spiral arm regions, 
through equation [21 by "str etching" the models to 
fit the observations fGG96. lMartmez-Garcfa et al 
2009ah : 



3.1. 



Azimuthal color gradient analysis 
with the GG96 method 



R„ 



Vrot 



I- age 



^'rad 

tani 



(2) 



The search and analysis of azimuthal color gra- 
dients are based on the four band, supergiant sen- 
sitive and reddening-fre^l photometric index 



Q{rJgi) = {r-J)- l^^ '^h g-i), (1) 
E{g - 1) 

where ^y- — 'Q- is the color excess term for a 
foreground screen. Q is a very good diag- 
nostic of star formation, since its value in- 
creases with the presence of blue and red su- 
pe rgiants. Details of the method can be found 



Rmeiai IS the mean radius of the studied region, 
Vrot is the galactic rotation velocitjQ from the lit- 
erature (corrected for inclination; see table [3]), 
d is the azimuthal distance from the shock (i.e., 
measured from d = 0), tage is the stellar model 
age at distance Urad is the radial velocity, 
and i is the spira l shock pit ch angle, assum- 
ing "steady state" (|Robertsl fl969) This equa- 
ti on can be easily derived from t he one presented 
in Martfnez-Garcfa et al. I (l2009al) . 



Gonzalez fc GrahamI d 19961 GG 96 hereafter). 



and Martmez-Garcfa et al. I ()2009al lbh. Briefly the 
procedure involves "unwrapping" the s piral arms 
by plotting them in a 9 vs. In i? map (e.g. . live et al 



19821 : lElmegreen et al.lll992r ). Under this geomet- 
ric transformation, logarithmic spirals appear as 
straight lines with slope = cot(-i), where i is the 
arm pitch angle. The search for color gradient 
candidates in the bar and arms can be done in 



*For a foreground screen, and for a mixtu re of dus t and 
stars as long as Ty < 2 l|Gonzalez fc Grahamj Il996l: 
iMartmez-Garci'a et all l2009al) . See also § 15.1.21 and fig- 
ure |48] 



^ For this research, we assume that i^i-ot ~ constant in the 
regions of interest. The mean value of i)rot for our sam- 
ple (after inclination correction, and excluding NGC 266 
due to its highly uncertain inclination angle) is 160 ± 30 
km s-l. The error in the velocity is comparable to ex- 
pected deviations from a flat rotation curve (e.g., se e the 
model rotation curves in iRomero-Gomez et al.ll2007D . We 
also try to avoid the inner parts of the disk, where a flat ro- 
tation curve should no longer be valid (the minimum value 
of -Rmcan/Rbar in figure [52] is ~ 0.4). 

■'For practical purposes, the d and tagc quantities introduced 
in equation [2] are really &d and 5tagc within the same az- 
imuthal region at a given radius (see also ^5.11 1. 

^ The term i^rad/ tan j, and equation [4] are only meaningful 
for spiral regions. 
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Table 3 
Galaxy parameters 



Name 


Type 


PA (deg) 


a (dog) 


Umax (km S ' ) 


Radial Velocity (km s ') 


Distance (Mpc) 


NGC 718 


SAB(s)a 


45 


29.4 ± 10.4 


58.5 ± 14.0 


1733 ± 10 


24.3 ± 2.1 


NGC 864 


SAB(rs)c 


20 


40.7 ± 3.1 


97.9 ± 4.1 


1560 ± 4 


22.0 ± 1.9 


NGC 4314 


SB(rs)a 


140" 


27.0 ± 5.3 


70.6 ± 5.0 


963 ± 26 


17.7 ± 2.1 


NGC 266 


SB(rs)ab 


95'' 


12.2 ± 16.5 


217.8 ± 7.5 


4661 ± 6 


64.8 ± 5.5 


NGC 986 


SB(rs)ab 


150 


40.7 ± 4.7 


43.3 ± 8.4 


2005 ± 7 


28.0 ± 2.4 


NGC 7496 


SB(s)b 


169. 7'^ 


24.2 ± 9.2 


65.0 ± 6.0 


1649 ± 6 


23.6 ± 2.0 


NGC 5383 


SB(rs)b 


85 


31. 7± 6.7 


142.6 ± 13.1 


2250 ± 4 


39.2 ± 3.4 


NGC 4593 


RSB(rs)b 


56'' 


42.2 ± 4.4 


161.4 ± 10.0 


2492 ± 6 


38.2 ± 3.3 


NGC 3059 


SB(rs)bc 


70.9= 


27.0 ± 8.4 


55.6 ± 5.5 


1260 ± 6 


15.0 ± 1.3 


NGC 7479 


SB(s)c 


25 


40.7 ± 3.1 


162.2 ± 7.1 


2378 ± 3 


34.4 ± 2.9 


NGC 3513 


SB(rs)c 


75 


37.4 ± 3.5 


38.6 ± 4.7 


1194 ± 7 


16.0 ± 1.4 



Note. — Col. (2). Hubble types from RC3. Col. (3). Position angles, mainly from RC3 and Hyperleda jPaturel et al.l 
120031) . Col. (4). Inclination angle, OL — COS {h/a)\ ajh is the isophotal diameter ratio derive d from the i?25 parameter in 
RC3. Col. (5). Maximum rotation velocity obtained from the HI data of IPaturel et ahl I I2003I) . uncorrected for inclination. 
Col. (6). Heliocentric radial velocity from RC3. Col. (7). Hubble distance obtained from the heliocentric radial velocity 
and the infall model of lMould et al.l 120001) . Ho = 71 ± 6 km s~' Mpc"^ 

° [Benedict et al.l 120021) . 

IPaturel eTall 120031) . 

" IPaturel et al.l 120031) : Lauberts A. & Valontijn E.A. (1989). 

''iVauglin et al.l |[T99gl). 



age processing compensates in pa rt for this ef- 



n„^-( f ?^Q;S^dt' - (0.hock + M)], (3) 



t 



R{t') 



A.^cotWln(||). (4) 

Equation [2] assumes nearly circular motion 
for the involved stellar regions. For the pat- 
tern speed determinations we take Wrad ^ 0, al- 
though this term beco mes important for tightly 
wound spirals (see also Grosbol fc Dottorl [2009'). 
Important deviations from circular orbits and 
velocity gradients are expected in the gas in 

the azimuthal direction, perpe n dicular to the 

bar touval fc Athanassoulall983 ). Martmez-Garci'a et al 
(j2009br ) have shown that, this notwithstanding, 
azimuthal color gradients across spiral arms can 
be detected, although assuming a circular motion 
dynamic model will result in a systematic trend 
to overestimate spiral pattern speeds at small 
radii, away from CR, in non-barred or weakly- 
barred galaxiesl3 Pixel averaging due to im- 



^^In these galaxies, gas flows to some extent along the arms 
after passing through a steady rotating spiral shock. The 
age gradient is narrower than in the absence of these non- 



fect ( Martfnez-Garcia et al.ll2009bl ). such that the 



measured pattern speed will be correct within the 
errors, but the systematic trend, whereby the dif- 
ference between measured and real pattern speeds 
DC , will still be observed. 

3.2. Bar extent 



According to lWozniak et al.l (|1995l ). the bar end 
is located after the radius of maximum ellipticity, 
and is marked by a change in the position an- 
gle (PA) of the isophotes. Isophotal PA, on the 
other hand, must remain roughly constant along 
the bar region, although spiral structure, rings sur- 
rounding the bar, and stellar bar ansae (mainly 



in ea rly-type galaxies, [Martinez- Valpuesta et al 



20071 ) may disturb the elliptical profiles. 

Although the bar extent may be underes- 
tima ted with the maximum ellip ticity crite- 
rion (jMichel-Dansac fc Wozniakl2006l) . this method 



circular motions, and hence an observer will infer a smaller 
difference between the orbital velocity of the stars and the 
pattern speed. Inside corotation, this will lead to an over- 
estimate of the pattern speed that increases as the radius 
decreases, and converges to zero at corotation, where the 
shock strength and the non-circular motions are minimal. 
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provides a good, homogeneous, standard for our 
purposes. We determined most bar lengths using 
the maximum elhpticity criterion (see table |4]). To 
this end, we masked bright stars and nearby ob- 
jects, and fitted elhpses to the bar's isophotes (in 
the NIR, mostly in the 2 ^m, deprojected frames) 
with the IR AF task ELLIPSE (|jedrzeiewskilll987t 
Buskolll99fi|) . We then generated plots of isophotal 
elhpt icity and PA versus radius (Iw ozniak et al 



that is absent in those with Mup 



lOMp 



1995: iMulchaev et alj Il997t iGadotti fc de Souza 



20061 ). The brightest pixel near the nuclear re- 



gion was given as an initial guess for the isophote 
center, while re-centering was allowed for outer 
isophotes. A second "mean isophote center" was 
computed with these re-centered isophotes; this 
second center was then used, without allowing for 
further re-centering. We also performed a visual 
estimate of the bar's extension. In figure [21 we 
compare the results of both measurements. 

Finally, we compare the bar length with the bar 
CR radius; the latter is calculated based on the 
pattern speeds derived from the color gradients in 
bar regions. 

4. Results 

In figures ordered by galaxy Hubble type, 
we show results for 42 regions with color gradi- 
ent candidates. Of these, 22 are located along the 
bar, and 20 are found in the spiral arms. For all 
objects, the direction of disk rotation was estab- 
lished assuming that spiral arms trail. Dust lane 
locations were inferred from the {g — J) color, and 
bar or spiral perturbations were traced in the NIR. 
We fit the data with the stellar population syn- 
thesis models of S. Chariot & G. Bruzual (2007, 
private communication); only models with solar 
metallicity are considered. The duration of the 
star formation burst is taken to be 2 x 10^ yr; a 
fraction of young stars of 2% by mass is mixed with 
a background of old stars 5 x 10^ yr old. We use a 
Salpeter initial mass function (IMF) with a lower 
mass limit Miowcr — O.IMq, and we try two dif- 
ferent IMF upper mass limits A/uppor: 10 and 100 
Mq (see figure [45|) . It is important to have in mind 
that real data may have Afuppor somewhere in be- 



tween these values. Models with Afu 
exhibit a sharp peak shortly after 



ppcr 



IOOM0 

= yr 



In real data this peak may be easily lost in low 
signal to noise regions or low resolution data, or 
confused with unnoticed artifacts or cosmic rays. 
Also, there appears to be an inverse correlation 
between gradient dctectability and emission 
(GG96. iMartmez-Garci'a et aDlioogal) . that could 
be explained if contamination from bright emis- 
sion lines produced in HII regions around the most 
massive stars masks the color gradients. We are 
able to fit a model with Muppcr = lOOAf© only in 
7 out of 42 regions (17%). 

For bar regions, we assume that stars age in the 
direction of rotation^ In spiral regions, we search 
for dust lanes upstream of star formation, and (or) 
establish the aging direction from the "match" be- 
tween the asymmetric profiles of the observations 
and the stellar models (see figure 05]). 

Each one of figures [3ll44l corresponds to an an- 
alyzed region. A deprojected mosaic of the galaxy 
(in the g-band, unless otherwise indicated), in 
logarithmic scale, is shown in the top left panel 
(a). The top right panel (6) displays the ob- 
served QirJgi) profile vs. azimuthal distance, d, 
in kpc {solid line and left y-axis); the observed 
(g — J) color vs. d (dotted line and second-from- 
right y-axis), with high values indicating dust lane 
locations; and the observed Kg surface bright- 
ness (in mag arcsec"^, lower values correspond to 
higher densities of old stars) vs. d [dashed line and 
rightmost y-axis). In the bottom left panel (c), 
a zoomed-in version of the Q{rJgi) vs. d profile 
(solid line) is compared with a stellar population 
model (dotted line) that has been "stretched" in 
iagc to fit the data. Model steUar age, tago, in units 
of 10^ yr, is shown in the upper x-axis. The verti- 
cal error bars show the size of photometric random 
errors, excluding the systematic calibration error; 
horizontal error bars represent possible deforma- 
tions of the Q(rJgi) profiles, due to density and 
metallicit y variations in the young a nd old popula- 
tions fsee iMartmez-Garcia et aL 20G9a) . We also 
compare the data with a mo del that cons iders the 
dissolution of stellar groups (jWielenlll977l ) after 50 



^*This is the aging direc tio n ex pected from ob- 
servat ions llZurita &: Perej l2008h and simula- 
tions jDobbs fc Pringldl2O10fK 



Ellipticity : 



minor axis 
major axis 
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Fig. 2. — Comparison between the bar extent obtained with the bar's isophotes (i?bar, see text), and the 
estimate by sight (i?visuai)- The dashed and dotted hnes indicate, respectively, the mean and the error of 
the ratio of both measurements, 1.02±0.04. The object numbers in the "New General Catalogue" (NGC, 
Dreverl[l888[ ). ordered by Hubble type, are shown in the horizontal axis. 
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Myr {dashed line^ see ^ 15.1. ip . Reduced values of 



2 / Qdata. Qmodcl 

X =^ . 



(5) 



were calculated for both models in the time inter- 
val -20 < iage(Myr) < 1000 The bottom right 
panel (d) exhibits isophote ellipticity vs. R^^^ in 
the upper plot, and isophotal PA vs. R^^^ in the 
lower one. Error bars were obtained from the EL- 
LIPSE task in IRAF. Hatched regions highlight 
the bar corotation region, as derived from the 
comparison between photometric data and stellar 
models. 

In table m we show the pattern speeds and reso- 
nance radii inferred from the comparison of stellar 
population synthesis models with observed color 
gradients. 

Remarks for each object: 

NGC 718 (Figures [3] - [5]). Regions A and B 
belong to the bar region. For both regions, the 
origin (tago — 0) of the stellar model that bests 
fits the observations is located in the leading side 
of the bar; the two regions yield a similar bar CR 
radius, within the errors. An inverse gradient is 
observed in spiral region NGC 718 C; this region 
was fit with Mapper = IOOMq. The arm and bar 
pattern speeds are similar in this object. 

NGC 864 (Figures i] - O. Although the po- 
sitions of color gradient candidates A and B, in 
the bar region, are quite different, relative to the 
bar Kg surface brightness and dust {g — J) pro- 
files, their analysis provides a similar CR position, 
within the errors. This barred galaxy has spiral 
arms with a "ragged" structure. Region C, lo- 
cated in the beginning of the eastern arm (left arm 
in the deprojected image), gives a corotation posi- 
tion near the end of the spiral arms, at ~ 76. 5 ±0.5 
arcsec. Region D, situated in an arm structure ap- 
parently decoupled from the main pattern, gives 
a corotation position similar to that of region C, 
within the errors. The two main arms of this ob- 
ject have different Q mean values, possibly owing 
to different levels of star formation activit\r3 



Although for some of the regions this time interval may in- 
clude structures not related to the color gradient (see, e.g., 
regions NGC 266 A & D, and NGC 3513 D; figures [T^ [TSl 
and 1421 respectively). 

This behavior was dubbed "Q effect" 
in lMartfnez-Garcfa et al.l ll2009af) . 



NGC 4314 (Figures [To] -E]). Two color gradi- 
ent candidates were found in the bar of this object. 
The one in region A is located in the leading side 
of the bar; the CR radius inferred from the com- 
parison between theoretical and observed QirJgi) 
profiles is near the location of maximum elliptic- 
ity. The color gradient candidate in region B is 
located in the trailing side of the bar; the inferred 
CR agrees, within the errors, with the result from 
region A. No color gradients were found in the 
arms. 

NGC 266 (Figures [H - [TS]). Unfortunately, 
the error due to the inclination angle, a, is 
higher than the a value itself (see table Al- 
though this produces very large errors in the com- 
puted r2p values, the errors in the resonance posi- 
tions are reason able, when equations Al and A3 
from Martfnez-Garcfa et al.l (2009a) are used. Re- 
gion A is the only one detected in the bar of this 
object, and its analysis yields a CR position close 
to the maximum of isophotal ellipticity; this region 
was fitted with a model with A/uppor = IOOMq. 
Regions B (presumably just before corotation, 
Mupper = IOM0) and D (presumably just after 
corotation. Mapper — IOOA/q) give similar reso- 
nance positions that also agree with the results 
of region A, within the errors. The resonance po- 
sition derived from region C is within 1 cr of the 
result of region A, and within 2 a of the position 
found from regions B and D. Regions B and C ap- 
pear to be associated with a "ring" feature, rather 
than with the spiral arms. 

NGC 986 (Figures [T6l-[T9l). This galaxy shows 
an important activity in the bar region when ob- 
served in the Q{rJgi) index. The bar's end is 
made evident more by the change in PA of the 
isophotes, than by the location of their maximum 
ellipticity. Region A harbors a color gradient can- 
didate in the bar whose analysis results in a CR 
radii close to the place where the isophotes change 
PA. Region B is located near the end of the bar; 



the fit of a stellar model (with 



IOOMq) 



to its observed Q{rJgi) profile gives a CR posi- 
tion that does not coincide with the bar's end. 
For regions C and D, in one of the arms, the coro- 
tation position coincides (within the errors) with 
the bar's end. 

NGC 7496 (Figures[20l-[24l). The northern spi- 
ral arm (lower arm in our deprojected image) dis- 
plays a compact elongated region of high surface 
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brightness, even in Kg] this arm is also more ex- 
tended in radius, when compared to the southern 
arm (upper arm in our deprojected image). Re- 
gions A and B are located in the bar and give sim- 
ilar CR radii, that encompass the location of max- 
imum ellipticity of the bar. Region C is probably 
located near the bar's CR, a fact that makes the 
determination of dynamic parameters uncertain. 
The observations, and the fit of the stellar model 
to regions D and E indicate inverse color gradients. 
Corotation is close to the ellipticity maximum of 
the bar; all the derived positions agree within 1 cr, 
except for the determination from region D. 

NGC 5383 (Figures - HH). This barred 
galaxy has very short spiral arms that do not reach 
the outer disk. The analysis of color gradient can- 
didates in regions A and B results in similar CR 
positions before the end of the bar. This may be 
due to strong non-circular motions that lead to an 
overestimate of ilp, and hence to an underestima- 
tion of the CR radius. 

NGC 4593 (Figures [H] - EH- Regions A 

(Mupper = IOOMq) and B (Muppor = lOMo), 

in the bar region, give a similar CR position near 
the maximum of ellipticity. Region C is located 
near the bar's end, and is probably associated 
with a "ring" feature. Region D, inverse color 
gradient in the spiral arm region, gives a similar 
corotation position in accordance with regions A 
and B. Region E ("direct" color gradient) yields a 
corotation position further away from the galaxy 
center than the one inferred from region D. The 
eastern and western sides of the bar (respectively, 
upper and lower sides in the deprojected frame) 
have different Q{rJgi) mean values. 

NGC 3059 (Figures [31] - [321) • This object is 
likely a double-bar system, and it is difficult to 
assess the bar's end from the ellipticity and PA 
vs. radius plots. However, the CR position derived 
from the bar region A lies close to the bar endpoint 
determined from visual inspection (^ 19.2 ±1.4 
arcsec). Region B in the arms, gives a corotation 
position at radii larger than the bar's end. 

NGC 7479 (Fi g ures p] - W)- According 



to iWozniak et al.l (|1995[ ). the HII regions lo- 
cated along the bar, near the dust lanes, show a 
"stretched" appearance probably due to strong gas 
flows that may trigger star formation. Regions A, 
B, (A/upper = IOMq), and C (Afuppcr - IOOMq) 
show candidate color gradients in the bar. The 



comparison of stellar models with observations 
renders CR positions close to the ellipticity max- 
imum (i.e., close to the bar's end). Region D, 
situated in one of the spiral arms and after the 
bar's corotation, features an inverse color gra- 
dient. The computed resonance positions agree, 
within the errors, for regions A through D. Region 
E lies near the bar's end; even though no color gra- 
dients are expected at this positionjlO the Q{rJgi) 
profile indicates that some star formation is taking 
place. 

NGC 3513 (Figures [31 -[HI). Regions B, C, E, 
and F exhibit several adjacent Q{rJgi) profiles, so 
it is hard to define a unique color gradient candi- 
date. Regions A and B belong to the bar, and their 
analysis locates CR just inside its end. Region C 
is located very close to the bar's end. Region D, in 
the spiral arm region, is presumably located inside 
corotation; the best fit of the models to the ob- 
served Q{rJgi) profiles suggests Muppor = lOOM© 
for this region, and yields a corotation radius close 
to the change in ellipticity and PA. In contradic- 
tion with the region D result, regions E and F in- 
dicate that the bar and spiral perturbations may 
rotate with different pattern speeds. Once again, 
Q{rJgi) has different mean values in both arms. 



At corotation the pattern and the rotating material have 
the same angular velocity, hence shocks should not occur 
(at least for non-barred spirals). 
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Table 4 

Observed and derived dynamic parameters. 



Galaxy and region 


Figure 


Location 


-Rmcan (arCSCc) 


flmoan (kpc) 


i^bar (arcscc) 




(km s ^ kpc ^ ) 


Ren (arcscc) 


RcR (kpc) 


NGC718 A 




3 




Bar 


20.7 ± 0.25 


2.4 ± 0.2 


22.4 ± 2.0 (ifs) 


2.6 ± 0.3 


37.9 ± 24.2 


26.7 ± 5.5 


3.1 ± 0.6 


NGC718 B 




4 




Bar 


21.2 ± 0.25 


2.5 ± 0.2 


22.4 ± 2.0 (Ks) 


2.6 ± 0.3 


41.4 ± 23.6 


24.4 ±3.6 


2.9 ± 0.4 


NGC718 C 




1 




Spiral 


41.4 ± 0.25 


4.9 ± 0.4 


22.4 ± 2.0 (k/) 


2.6 ± 0.3 


33.6 ± 13.2 


30.1 ± 3.0 


3.5 ± 0.4 


NGC864 A 




6 




Bar 


26.0 ± 0.25 


2.8 ± 0.2 


36.1 ± 3.5 {K,) 


3.8 ± 0.5 


38.6 ± 6.5 


36.5 ± 4.7 


3.9 ± 0.5 


NGC864 B 




7 




Bar 


27.7 ± 0.25 


3.0 ± 0.3 


36.1 ± 3.5 (K^) 


3.8 ± 0.5 


44.1 ± 6.6 


31.9 ± 3.4 


3.4 ± 0.4 


NGC864 C 




1 




Spiral 


40.4 ± 0.25 


4.3 ± 0.4 


36.1 ± 3.5 (K^) 


3.8 ± 0.5 


18.6 ± 4.1 


75.5 ± 13.5 


8.1 ± 1.4 


NGC864 D 




9 




Spiral 


61.1 ± 0.25 


6.5 ± 0.6 


36.1 ± 3.5 {Ks) 


3.8 ± 0.5 


17.0 ± 3.1 


82.6 ± 11.0 


8.8 ± 1.2 


NGC4314 A 




10^ 


Bar 


44.6 ± 0.25 


3.8 ± 0.5 


70.3 ± 17.0 {Ks) 


6.0 ± 1.6 


30.6 ± 10.4 


59.3 ± 11.0 


5.1 ± 0.9 


NGC4314 B 


m 


Bar 


30.4 ± 0.25 


2.6 ± 0.3 


70.3 ± 17.0 {Ks) 


6.0 ± 1.6 


42.4 ± 15.5 


42.8 ± 8.8 


3.7 ± 0.8 


NGC266 A 


ll2| 


Bar 


16.0 ± 0.14 


5.0 ± 0.4 


17.5 ± 1.7 (K') 


5.5 ± 0.7 


193.1 ± 452.9 


17.0 ± 2.4 


5.3 ± 0.8 


NGC266 B 




SpiraF 


19.6 ± 0.14 


6.2 ± 0.5 


17.5 ± 1.7 (K') 


5.5 ± 0.7 


158.3 ± 372.7 


20.7 ±3.3 


6.5 ± 1.0 


NGC266 C 


M 


SpiraF 


15.1 ± 0.14 


4.7 ± 0.4 


17.5 ± 1.7 (K') 


5.5 ± 0.7 


213.0 ± 485.6 


15.4 ± 1.8 


4.8 ± 0.6 


NGC266 D 


[m 


Spiral 


22.9 ± 0.14 


7.2 ± 0.6 


17.5 ± 1.7 Ik') 


5.5 ± 0.7 


158.6 ± 318.1 


20.7 ± 2.7 


6.5 ± 0.8 


NGC986 A 


lisi 


Bar 


26.9 ± 0.15 


3.7 ± 0.3 


50.5 ± 4.6 {Ks)"" 


6.8 ± 0.9 


11.2 ± 4.3 


43.9 ± 9.2 


6.0 ± 1.3 


NGC986 C 


EI] 


Spiral 


52.4 ± 0.15 


7.1 ± 0.6 


50.5 ± 4.6 (ifs)" 


6.8 ± 0.9 


20.5 ± 2.5 


23.8 ± 2.8 


3.2 ± 0.4 


NGC986 D 


EI] 


Spiral 


68.9 ± 0.15 


9.4 ± 0.8 


50.5 ± 4.6 (ifs)" 


6.8 ± 0.9 


11.4 ± 1.9 


42.8 ± 4.0 


5.8 ± 0.5 


NGC986 E 


EI] 


Spiral 


83.7 ± 0.15 


11.4 ± 1.0 


50.5 ± 4.6 (ifs)" 


6.8 ± 0.9 


10.1 ± 1.6 


48.2 ± 4.7 


6.5 ± 0.6 


NGC7496 A 




201 


Bar 


21.5 ± 0.15 


2.5 ± 0.2 


38.0 ± 3.5 {Ks) 


4.3 ± 0.5 


39.6 ± 32.5 


35.0 ± 14.4 


4.0 ± 1.6 


NGC7496 B 


[21] 


Bar 


23.0 ± 0.15 


2.6 ± 0.2 


38.0 ± 3.5 {Ks) 


4.3 ± 0.5 


34.6 ± 30.2 


40.1 ± 18.9 


4.6 ± 2.2 


NGC7496 C 


[22] 


Bar 


33.0 ± 0.15 


3.8 ± 0.3 


38.0 ± 3.5 {Ks) 


4.3 ± 0.5 


30.1 ± 20.2 


46.0 ± 11.2 


5.3 ± 1.3 


NGC7496 D 


[23] 


Spiral 


44.1 ± 0.15 


5.0 ± 0.4 


38.0 ± 3.5 {Ks) 


4.3 ± 0.5 


51.4 ± 13.8 


27.0 ± 4.2 


3.1 ± 0.5 


NGC7496 E 


[21 


Spiral 


68.4 ± 0.15 


7.8 ± 0.7 


38.0 ± 3.5 {Ks) 


4.3 ± 0.5 


31.1 ± 9.3 


44.6 ± 5.9 


5.1 ± 0.7 


NGC5383 A 
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Bar 


20.1 ± 0.25 


3.8 ± 0.3 


54.2 ± 5.0 {Ks) 


10.3 ± 1.3 


65.2 ± 19 


21.9 ± 2.4 


4.2 ± 0.5 


NGC5383 B 


[26] 


Bar 


27.3 ± 0.25 


5.2 ± 0.5 


54.2 ± 5.0 {Ks) 


10.3 ± 1.3 


48.8 ± 13.6 


29.2 ± 3.0 


5.6 ± 0.6 


NGC4593 A 


[27] 


Bar 


34.8 ± 0.14 


6.4 ± 0.6 


48.0 ± 4.4 {K) 


8.9 ± 1.1 


21.8 ± 6.4 


59.5 ± 12.9 


11.0 ± 2.4 


NGC4593 B 


[28] 


Bar 


36.7 ±0.14 


6.8 ± 0.6 


48.0 ± 4.4 {K) 


8.9 ± 1.1 


26.8 ± 5.6 


48.5 ± 6.6 


9.0 ± 1.2 


NGC4593 C 


[29] 


Spiral^ 


56.1 ± 0.14 


10.4 ± 0.9 


48.0 ± 4.4 {K) 


8.9 ± 1.1 


10.3 ± 4.1 


125.8 ± 41.4 


23.3 ± 7.7 


NGC4593 D 


[30] 


Spiral 


72.6 ± 0.14 


13.4 ± 1.2 


48.0 ± 4.4 {K) 


8.9 ± 1.1 


27.0 ± 2.8 


48.1 ± 3.0 


8.9 ± 0.6 


NGC4593 E 


[31] 


Spiral 


64.1 ± 0.14 


11.9 ± 1.0 


48.0 ± 4.4 {K) 


8.9 ± 1.1 


11.1 ± 3.3 


117.0 ± 25.7 


21.7 ± 4.8 


NGC3059 A 


[32] 


Bar 


12.4 ± 0.14 


0.9 ± 0.1 


19.2 ± 1.4 {H)" 


1.4 ± 0.2 


121.2 ± 49.8 


13.9 ± 1.6 


1.0 ± 0.1 


NGC3059 B 


[33] 


Spiral 


26.5 ± 0.14 


1.9 ± 0.2 


19.2 ± 1.4 {H)'' 


1.4 ± 0.2 


58.6 ± 23.6 


28.7 ± 3.1 


2.1 ± 0.2 


NGC7479 A 






Bar 


25.0 ± 0.25 


4.2 ± 0.4 


54.2 ± 5.0 {Ks) 


9.0 ± 1.1 


29.5 ± 8.1 


50.5 ± 11.3 


8.4 ± 1.9 


NGC7479 B 


[35] 


Bar 


47.6 ± 0.25 


7.9 ±0.7 


54.2 ± 5.0 {Ks) 


9.0 ± 1.1 


19.8 ± 4.3 


75.2 ± 12.3 


12.5 ± 2.1 


NGC7479 


[36] 


Bar 


41.4 ± 0.25 


6.9 ± 0.6 


54.2 ± 5.0 {Ks) 


9.0 ± 1.1 


21.9 ± 4.8 


68.2 ± 11.6 


11.4 ± 1.9 


NGC7479 D 


[37] 


Spiral 


80.3 ± 0.25 


13.4 ± 1.1 


54.2 ± 5.0 {Ks) 


9.0 ± 1.1 


27.5 ± 2.2 


54.2 ± 3.3 


9.0 ± 0.6 


NGC7479 E 


[38] 


Spiral 


60.1 ± 0.25 


10.0 ± 0.8 


54.2 ± 5.0 {Ks) 


9.0 ± 1.1 


14.4 ± 3.4 


103.7 ± 19.1 


17.3 ± 3.2 


NGC3513 A 
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Bar 


14.8 ± 0.14 


1.1 ± 0.1 


24.6 ± 2.5 {H) 


1.9 ± 0.3 


45.0 ± 10.1 


18.2 ± 2.2 


1.4 ± 0.2 


NGC3513 B 


\M 


Bar 


17.0 ± 0.14 


1.3 ± 0.1 


24.6 ± 2.5 (H) 


1.9 ± 0.3 


39.5 ± 8.9 


20.7 ± 2.5 


1.6 ± 0.2 


NGC3513 C 


\M\ 


Bar 


24.1 ± 0.14 


1.9 ± 0.2 


24.6 ± 2.5 (H) 


1.9 ± 0.3 


28.2 ± 6.1 


29.0 ± 3.4 


2.3 ± 0.3 


NGC3513 D 


\M 


Spiral 


20.8 ± 0.14 


1.6 ± 0.1 


24.6 ± 2.5 {H) 


1.9 ± 0.3 


22.9 ± 7.3 


35.8 ± 7.5 


2.8 ± 0.6 


NGC3513 E 


\M 


Spiral 


43.0 ± 0.14 


3.3 ± 0.3 


24.6 ± 2.5 (H) 


1.9 ± 0.3 


9.9 ± 3.7 


82.7 ± 19.8 


6.4 ± 1.5 


NGC3513 F 


[44] 


Spiral 


35.8 ± 0.14 


2.8 ± 0.2 


24.6 ± 2.5 (H) 


1.9 ± 0.3 


14.7 ± 4.2 


55.6 ± 9.8 


4.3 ± 0.8 



Note. — Col. (4) and (5). Mean radius of the studied regions, in arcsec and kpc, respectively. Col. (6) and (7). Radius of the bar, in arcsec and kpc, respectively, and 
bandpass used to determine it. Col. (8) Pattern speed. Col. (9) and (10). Corotation radius, in arcsec and kpc, respectively. 

^Obtained from the change of PA with radius of the bar isophotes. 

^Obtained from visual estimate. 

^These regions seem to be associated with rings, rather than with spirals. 



5. Discussion 

5.1. Characteristics of observed color gra- 
dients 



In a "standard color gradient picture" , assum- 
ing circular motions, one would expect an az- 
imuthal sequence starting with compression of gas, 
dust lanes, star formation onset and stellar drift, 
which would lead to color gradients. When one 
adds an extended period of star formation, disper- 
sion velocities and post shock velocities, the pre- 
dict ed color profiles niay be come rather broad (see, 
e.g- lYuan fc GrosbollflQSl . 

Although shocks and/or dust lanes are located 
mainly upstream relative to the spiral poten- 
tial minimum, in a more detailed scenario, for 
certain models and relatively short radial in- 
tervals, they can be found downstream (in the 
gas st ream direction) the sp iral potential mini- 
mum (|Gittins fc Clarkell2004 ). This is due to the 
fact that pitch angles of the arms are expected 
to follow ip > Id > isF, where ip is the pitch 
angle of the potential, ip, is the pitch angle of the 
dust lane (the shock) , and isp is the pitch angle of 
star formation. According to this, the azimuthal 
relative location of potential and dust lanes may 
change within the same galaxy and be radially 
dependent. In principle, the local potential min- 
imum could be determined from the the K-hand 
spiral arm0 Nevertheless, gravity is a long-range 
force, and all the non-axisymmetric contributions 
must be taken into account for local potential min- 
imum determinations. Thus, the local potential 
does not necessarily coinci de with the l o cal den 
sity of old stars (see, e.g., 'Zhang' 'l996'; 'Berma: 
20031 : [Cittins fc Clarke 2004; Zhang & Buta 20o' 



Buta fc Zhanglfeooit ). and the expected locations 



of the dust lanes must follow the potential, rather 
than the if-band intensity. 

For some of the bar regions presented in fig- 
ures [3]l44l the tago = location of the gradient 
seems to be located before the main dust lane 
peak. This can be attributed to the fact that star 
formation in bars is supposed to begin at the dust 



^*As already mentioned in §[3l although young stellar popu- 
lations may ac count for only 3% of the global i^-band flux 
iRhoadslllQQSl ). they may contribute up to a third of the 2 
lj,m emission in local features. 



"spurs" of the bar itself (see § 11.11) , which are lo- 
cated in the "trailing" side and upstream the main 
dust lanes. 

For this investigation the d — distance (i.e., 
the assumed shock location) is chosen where the 
maximum in the g — J color is observed. This 
is the case, even for regions where double peaks 
or dust lanes are seen (e.g., region NGC 864 B, 
figure [7]). For most of the fits we do not have 
^agc = for d = precisely (e.g., region NGC 
7496 A in figure [20|. This is due to the fact that 
only the width of the model profile in Q{rJgi) is 
fitted to the observations by stretching. Theoret- 
ically, for spiral regions tage = should coincide 
or be located after d = 0. In real galaxies, it is 
hard to pinpoint where the onset of star forma- 
tion really occurs. One important reason is the 
fact that a diffuse cloud of neutral gas has first to 
become a dense cloud, then a molecular cloud, and 
finally a self-gravitating cloud to achieve star for- 
mat ion; this process may take ~ 10^ years (see, 
e.g.. iTamburro et al. I I2OO8I: lEgusa et alj f l2009h . 



5.1.1. Downstream decline of the gradients 



Martmez-Garcia et al. ( 2009a[ ) already noticed 
that for some spiral regions there is a "downstream 
decline" (or "downstream fall") of the gradients. 
In such regions, the observed Q(rJgi) profile de- 
clines below the model (or the value Q{rJgi) ^ 
1.57, for solar metallicity, in figure I45|) . appar- 
ently returning to the "old background popula- 
tion" level, or lower, by i « 5 x 10^ years, much 
faster than the theoretical expectations. This is 
the case, for example, of regions NGC 864 C, NGC 
7496 B & C, NGC 7479 B, NGC 3513 A & C 
(figures El EB Ea Ea [391 and [HI respectively), 
for which the fit between Q{rJgi) data and dot- 
ted line model is not good, except for the range 
< tage < 50 Myr 
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Martinez-Garcia et al. ( 2009al ) hypothesized 



that the decline of the observed Q{rJgi) profiles 
below the models (that assume pure circular or- 
bits) might be caused by s tellar non-circular mo- 
tions in the data. However. iMartmez-Garcia et al. 
( 2009bl ) found that non-circular motions cannot 



^^Incidentally, Bruzual & Chariot models previously to 1997 
agreed much better with this rapid decline of Q; see GG96, 
their figure 19. 
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explain the discrepancy. 

A solution to this problem may be provided by 
the dissolution of stellar groups scenario proposed 
bv lWieleii ( 1977 . see figure According to this 
author, the diffusion of stellar orbits can enhance 
the dissolution of young stellar groups by increas- 
ing their internal velocity dispersion. To explain 
the observed in crease in the v elocity dispersion of 
stars with age, IWielenI ()l977l) proposes the exis- 
tence of local fluctuations of the gravitational field 
with a rather stochastic behavior. This irregular 
field has the effect of creating a diffusion process 
in the velocity space of stars. The dissolution of 
stellar groups proceeds in two phases. During the 
first phase, the internal velocity dispersion, tTint, 
causes the group to expand, until its diameter, 
Dsg(iagc) is larger than the distance over which si- 
multaneous perturbations (for stars close together 
in space) are significantly correlated, i.e., the co- 
herence length, Lco0 During this first phase, the 
center of mass of the group is affected by diffusion 
mechanisms, while the members of the group only 
suffer small tidal effects, thus aint ^ constant. The 
duration of the first phase is given by: 



^'phascl 



int ; 



(6) 



For typical valu es of gint =10 km s ^, 
- I?,g(0)=500 pc (|Wielenlll977l) . we have 



for tphasei in Myr, Leo - -Dsg(O) in pc, and crint in 
km s~^. 
and Leo 
that: 

iphasci ^ 50Myr. 

During the second phase of the dissolution (af- 
ter tphasci years), the perturbations over each star 
member of the group are not longer correlated, 
and the group dissolves with time because of the 
diffusion mechanism. During this phase the veloc- 
ity dispersion increases with time, i.e., a — <y{t), 



where t = t^ 



phasel • 



In order to evaluate the relevance of the dis- 



solution of stellar groups scenario (jWielenI 119771 ) 
for the color gradient picture, we have built a 



^"This length depends on the mechanism that causes 
the irregular gravitational field, with the consequent 
disk heating. Possible mechanisms may be provided 
by g iant molecular clouds, spiral arms llLacevI 1 1984 



1991 ) , small-scale dark matter clumps llBerezinskv et al] 



20031: iBarranco fc Bernall l201lD . or massive dark clus- 
ters dSanchez-Salceddll99Sr " in dark matter halos. 



"toy model" for Q{rJgi), consisting of two phases. 
During the first phase, while tage < 50 Myr, the 
fractions by mass of young, /3i = 2%, and old, 
— 98%, stars are kept constant. During the 
second phase, we assume (1) that the fraction of 
young stars changes as /3i oc l/^ago: (2) that 
the surface density of y oung stars be haves in the 
same way (see eq. 10 in lWielen|[l977l) . 

The models produced with this approximation 
are displayed in figure S?! with a continuous line 
for models with the IMF Afuppcr — IOMq, and 
a dotted line for those w ith Afupper = 100Mf7) 
(see al so figures 8 and 9 in iMartfnez-Garcfa et al 
2009al ). According to figure H?! the dissolution 



of stellar groups, as a consequence of disk heat- 
ing, may provide an efficient mechanism to ex- 
plain the observed "downstr eam decline" of the 



gradien ts encountered in Martfnez-Garcfa et al 



(|2009al ). and in this investigation. It is important 
to mention that the dissolution effect operates also 
in the presence of non-circular motions. 



Other dissolution 



scenarios: 



Another effective disruption mechanism for star 
clusters comes from stellar winds and supernovae 
explosions that remove gas (and dust) on short 
timescales. These perturb the potential and cause 
young star clusters to become unbound (see, e.g., 
Bastian fc Goodwin|[200l iGieles fc Bastianll2008l: 
Lamers et al.l 2005 ). This "infant mortality" oc- 
curs during the first 10 Myr and affects the most 
luminous clusters (those with the most rapid color 
evolution). Nevertheless, stars escape with the ini- 
tial velocity dispersion of the cluster and are phys- 
ically a ssociated with it for 10-40 Myr after gas re- 
moval ( Bastian fc Goodwin l2006l) . The color gra- 
dients presented in this investigation are probably 
formed by clusters that survive "infant mortality" 
at least for 50 Myr. 

5.1.2. Color gradients in dusty environments 



One more issue in the "standard color gradi- 
ent picture" relates to the Q(rJgi) value expected 
for the old background population at iagc < 0. 
Also, Q{rJgi) (or any color) should return to the 
background value on both sides of the bar, at 
the same spatial offsets. In some of the observed 
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color gradient candidates in this investigation, the 
Q{rJgi) index does not agree well with the mod- 
els for tage < (see, e.g., regions NGC 718 C, 
NGC 986 A, NGC 7479 A; figures [5l [H and El 
respectively) . 

Models with pure old background population 
have an approximately constant Q{rJgi) value 1^ 
regardless of surface brightness. However, the ex- 
act value depends on the av erage age and metallic 



ity of the local region (see iMartmez-Garcia et al 



2009al ). that cannot, especially metallicity, be de- 



termined unambiguously just from photometric 
data. A complementary possible explanation of 
why, for some regions, the QirJgi) index does not 
decline to the background value for iagc < is 



that the adopted method (see § 13. 1[) involves aver- 
aging over radial annuli to increase the S/N ratio. 
Albeit we do so carefully, we may be combining 
zones with somewhat different background levels. 

Yet another recurring concern is whether dust 
can mimic the Q profile of a star formation burst, 
since Q{rJgi) is reddening-insensitive for a fore- 
ground screen, but not exactly so for a mix- 
ture of dust and stars ( Gonzalez &: Graham 19961 : 



iMartmez-Garcia et al ]|2009ah . The main expected 
effect would be a higher Q{rJgi) value (see iig- 
ure 148)1 . 

19991) 



Assuming ry ~ 1.0 (e.g., IXilouris et al 



for a nearly face-on disk, given the inclination an- 
gles of the galaxy sample, our observations cover 
on average the range < ry < 2.0. For regions 
where tagc < 0, though, thick dust lanes may be 
present with ry 2.0. 

From the lCharlot fc Falll (|2000l see also Bruzual 
& Chariot 2003) dust model, we get the linear re- 
lation Q{rJgi) = 0.02tv -|- 1.51 (where 1.51 is 
the background Q{rJgi) value; see figure [48] for 
iagc ^ 0)o Thus, to get a Q{rJgi) profile rem- 
iniscent of a density wave induced color gradient 
in the absence of young stars, with a peak value 
^ 1.66, the dust must have ry ~ 7.0. But also, for 
this scenario to take place, ry must increase for 



Q{rjgi) only stabilizes after a couple of Gyr. Although we 
are using a background population 5 Gyr old, in actuality 
stars take at most a few hundred Myr to go from one arm 
to the next, so that the background will have a contribution 
from several previous bursts between leS and le9 yr old, 
for which Q is higher and still slowly declining. 

The mean photometric error in Q (excluding the zero point 
error) is 4(TQ(^jgi) < 0.06 mag. 



Q(rJgi) to increase, and viceversa, i.e., the profile 
in a reddening sensitive index, like {g — J), must 
follow the Q{rJgi) azimuthal distribution. 

For iagc > 0, dust becomes progressively less 
important, because star formation processes (e.g., 
UV radiation) sweep away and destroy available 
material. Indeed, concentrations of dust that are 
higher downstream the shock than at the shock are 
not observed in most real spiral arms. Fittingly, 
all of the observed color gradients indicate an in- 
verse correlation between Q[rJgi) and {g — J), so 
that dust cannot be mimicking a star formation 
burst El Although the situation may be different 
for bar regions, we do not observe a correlation be- 
tween the dust and Q{rJgi) profiles in any of the 
bar gradient candidates analyzed in this investiga- 
tion. 

Finally, because of the connection between the 
gradients (i.e., star formation) and disk dynamics 
(i.e., the orbital resonance positions; see § 15. 2p . 
we argue that dust features are not affecting the 
analyzed gradients in any important way. 

5.2. Connection with dynamics 



In figure I49( the pattern speed obtained from 
each region in NGC 7479 is plotted vs. the mean 
radius of the region. Notice the tendency to de- 
rive slower Vl!^'" from bar regions (A, B, C) at 
larger radius. Likewise, figure [50] shows the fi^^'' 
obtained for NGC 3513, vs. the radius of each re- 
gion. As in plot |49] (NGC 7479), a trend whereby 
rip'^'' cx i s observed for the bar regions (A, 
B, and C). In iMartfnez-Garcia et al.l (|2009b[ ). we 
found that this effect is observed in non-barred 
or weakly barred galaxies, if non-circular motions 
are present but color gradient data are analyzed 
assuming stars have purely circular orbits. In the 
case of barred galaxies, non-circular motions in the 
bar cause fi^^'' to be overestimated inside the bar 



CR radius; once more, the size of the discrepancy 
between the real and the measured pattern speed 
increases as the radius decreases!^ This signature 



■^•^ In this regard, Gonzalez ct al. (1996) found that the gradi- 
ent in M99 also follows the models when mapped in colors 
sensitive to dust, such as {g — i), {g — K), and {g — J). This 
would not be the case if the profiles were contaminated by 
dust. 

■^■^ Again, to explain a "thinner" gradient the observer has to 
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does not mean that fi^'^'' is changing with radius, 
in either kind of galaxvl^ 

Given that many of the computed resonance 
positions seem to be in accordance with theoret- 
ical predictions (e.g., bars ending near their CR 
radius), most of the analyzed color gradient can- 
didates may indeed have a relation with the dy- 
namics of the disk. Figure [51] shows a plot of 
bar extent, i?bar (see § 13. 2p . vs. bar CR radius, 
i?CR) as inferred from the comparison between 
stellar population models, and broad-band opti- 
cal and NIR observations. Red open triangles de- 
note bar regions, whereas black solid circles indi- 
cate color gradient candidates in the arms. The 
plot is divided in three zones, i.e., the "slow", 
"fast" , and "forbidden" bar areas; the latter corre- 
spond s to the "super-fast" bars of Buta fc Zhand 
(|2009l ). The dotted line, Rcr = Rhar, and the 
dashed line, Rcr — 1.4i?bar, enclose the param- 
eter space where most of the points should be 
if bars end near CR ( Aguerri et al.l [200 3). If 



the spiral pattern speed, n^P'^''^^ is similar to 
the bar's pattern speed, ft^'^'^, then spiral re- 
gion points should fall in this area, too. Un- 
fortunately, with the GG96 method it is difficult 
to distinguish between "super-fast" bars and the 
expected overestir nation of ft^'^'^, owing to non - 



circular motions ( Martmez-Garci'a et al. I l2009bl) 
All the points below the Rcr = Rhai (dotted) line 
are more likely due to the latter effect. In order 
to dis criminate between these two po ssibilities, we 
follow iMartmez-Garci'a et al.l (|2009bh and define: 



data 



(7) 



Here, Odata is the pattern speed value obtained 
with the GG96 method (either for the bar or the 
spiral), and fip is the pattern speed of the (bar 
of spiral) perturbation derived from the rotation 
curve once a resonance position is fixed. 

For barred spirals, in the case where bar and 
arms share the same pattern speed (and thus the 
same CR radius), we have: 



Vrot 



1.2i?ba 



(8) 



invoke a smaller difference between the orbital and the pat- 
tern speeds; since orbital speed increases inward, so must 
the measured S7p^''. 
■^^If it were, it would not be a pattern speed! 



where we adopt ^ 1.2, the expected ra- 

tio for spirals with well defined bars (se e, e.g., 
lAthana ssoula 1992; Elmcgrccn 1996; Buta fc Zhang 
I2OO9I ). In figure[52l we show Sil^ vs. i?mcan/l-2-Rbar; 
this ratio is ~ 1 for regions near CR. Once again, 
bar regions are shown with open red triangles to 
distinguish them from spiral regions (solid black 
circles). For the bar regions, there is a system- 
atic trend, whereby fl}^^'^ is overestimated at small 
radii, the magnitude of the bias is inversely cor- 
related with radius, and the points converge to 
6flp = as they approach the corotation zone. 
This is the expected effect when non-circular mo- 
tions are present, but color gradients are inter- 
preted with a dyna mic model that considers purel y 
circular motions ( Martmez-Garcia et al. l2009bf ). 
However, as a consequence of data processing that 
was already discerned by these authors, regions in 
the outer half of the bar yield fi^^'' values with 
less than 50% error. The detectable trend, on the 
other hand, confirms the link of the gradients to 
disk dynamics. 

In figure [53l we plot Sflp for the spiral arm re- 
gions vs. region radius, in the same units as fig- 
ure [52I we highlight those with "inverse" color 
gradients (i.e., the sense of rotation is opposite 
to the observed aging of stars) with open red tri- 
angles. Most of these regions lie over the dotted 
line, where the pattern speed obtained from the 
color gradient candidates, fidata, equals the pat- 
tern speed derived from a flat rotation curve, if 
the bar CR is located at 1.2 i?bar- Conversely, 
regions where stars age in the direction of disk 
rotation (solid black circles) sit below the dotted 
line. If such "inverse" color gradients occur in ob- 
jects where 51^'^'' « fl^"^^, then this plot may in- 
dicate the presence of non-circular motions close 
to corotation, that cause the overestimation of the 
pattern speed with our method. 

On the other hand, the existence of gradients 
where stellar aging follows disk rotation, at radii 
beyond the bar CR radius, may be interpreted as 
7^ OpP"''*^ i.e., decoupled pattern speeds for 
the bar and the spiral; the CR radius of the spiral 
pattern would be larger than that of the barrel 



26 Two of these regions (NGC 266 B & C) seem, at least 
visually, more associated with rings rather than with spiral 
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5.3. The origin of spiral arms. 



The origin of spiral arms may be differ- 
ent in barred galaxies than in non-barred or 
weakly-barred spirals. An alternative theory to 
density waves propounds that manifold- driven 
chaotic orbits are the foundation of spirals and 
rin gs (nuclear rings excluded) in barred galax- 
ies (|Patsis"2006'; 'Romero- Gomez et al."2006'. 2007t 



Athaiiassoula ct al. 2009allbl. 12010: Voglis et al 
2006ai bl: lTsoutsis et al.ll2008l 



2009; 



Harsoula fc KalapotharaEog 



gin of the arms and/or their subsequent evolution 
may follow different paths from the ones proposed 
by the "manifold" theory. 

Concerning a different but related aspect, the 
observed trends for color gradients with radius in 
figure [53] suggest that non-circular motions are im- 
portant near the bar corotation in barred galaxies, 
and that their effect on the determination of the 
spiral pattern speed are significantly stronger for 
inverse color gradients, i.e., when ilp'^'' Qsp"-ai|27| 
The [p resence of non-circular motions is not sur- 



onnnkr'^ ■ TTT n ' + i ' /'onnci prTsftg, given the existence of the bar. If, however 

2009 1. In the view of Romero-G omez et al. (20061 ^, • r rr^ ■ i i i • i -^i 



20071 ) and lAthanassoula et all (2009a. b. 2010/), 



invariant manifolds are "tubes" that guide or- 
bits; they are associated with unstable Lagrangian 
points at corotation. Material is trapped in 
the manifolds during disk evolution, and cir- 
culat es across the disk, ind ucing radial mix- 
ing ( Athanassoula et al. 2010l) . Another inter- 



pretation of the "invariant manifold theory" con- 
siders only apsidal (apocentric or pericentric) sec- 
tions of unstable manifold s ((Voglis et al. 2006a.bl: 
Tsoutsis et al1l2008l liopgl ). In this latter version. 



there is no need for c onstantly supplying material 
inside the manifolds (Eft hvmiopoulos. C.l 2010[ ) . 

One important prediction of the "manifold the- 
ories" (both views) is that the spiral arms and the 
bar must have the same pattern speed. Accord- 
ing to the age gradients found in the spiral arms 
of our sample, the objects NGC 718, NGC 266, 
NGC 986, NGC 7496, NGC 4593, and NGC 7479 
are the most likely to have rjj^^'' « rj^?'"^^'. Con- 
versely, NGC 864, NGC 3059, and NGC 3513 
appear to have different spiral and bar pattern 
speeds. No conclusion can be drawn for NGC 4314 
and NGC 5383, since we were able to detect gra- 
dients only in the bars of these objects. Quite 
interestingly, with the exception of NGC 7479, 
all objects in our sample with VH^'^ ss ri^P"'''' 
are early Hubble types, while the objects with 

fibar ^ f^spiral j^^^g Hubblc typCS. 

Both the density wave and the "manifold" the- 
ories predict a spiral potential that can produce 
shocks in the circulating gas. If star formation is 
triggered by these shocks, it would be difHcult to 
discriminate between the two theories on the ba- 
sis of the presence of color gradients alone. But, 
at least in the three objects in our sample with 
uncoupled bar and spiral pattern speeds, the ori- 



the size of (5f2p in barred galaxies correlates with 
shock strength, the amplitudes of coupled spiral 
arms should also be stronger near CR (see, e.g. 



NGC 1566 spiral Fourier amplitudes in lSalo et al 
201fl[ ). Stron g shocks are neither expected nor 



observed (e.g.. lMartfnez-Garcfa et al. 2009b) near 
the pattern CR radius of non-barred or weakly- 
barred spiral galaxies. 



After a re-analysis of the result in iButa et al 



(|2009[ ) who find only a weak indication that some 



stron; 



(201 
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bars may drive strong spirals, iSalo et al 
conclude, contrariwise, that in a statis- 
tical sense spiral density waves may indeed be 
driven by bars. More comparisons w ith observa- 
tions (see, e.g., ICrouchv et al. 20101 ) are needed 
to test b oth the "manifold" and th e density wave 
theories ( Athanassoula et al. 2010[ ). 

To confirm the link between the gradients and 
the disk dynamics, it is also important to compare 
the location of the spiral " end point s 'r^l \v ith the 
orbital resonance positions. Schward ( 19851 ) shows 
that bar-induced gas spiral arms may reach be- 
yond the outer Lindblad resonance (OLR) of the 
bar. In figure [M] we plot the OLR of the bar, 
-^OLR' assuming 



^^For inverse color gradients, 5f!p = 1.55(1.2il|3g^j./_Rinean) — 
1.07 (pearson correlation coefficient, p = 0.51); for gradi- 
ents with aging in the rotation direction (i.e., decoupled 
pattern speeds), (5r2p = 0.89(1.2_Rbar/fimcan) — 1-19 (p = 
0.73). 

■^^ Strong bars or spirals may be defined by comparing 
the non-axisymmetric gravitational perturbation (that in- 
duces a ta ngential force) to the mean axisymmetric radial 
force (e.g.. lBlock et al.ll200j) . 

■^^Or the maximum radial extent of the arms, since spirals 
may fall ba ck towards smaller radii l|Athanassoula et al.l 
l2OO9bl . [2O10l) . 
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R 



bar 

OLR ' 



1.2i?ba 



V2' 



(9) 



vs. the spiral extent, i?™, estimated by eye in the 
NIR, for ah objectf[f^ in our sample (arms with 
decoupled pattern speeds are shown as red open 
triangles). R^™^ values are listed in table O 

The dotted line indicates the identity, -R^,™ = 
i?OLRj while the dashed line is the OLS (ordi- 
nary least squares) bisector On average, for our 
whole barred spiral galaxy sample, the radial ex- 
tent of the spiral arms coincides with the OLR of 
the bar, regardless of the pattern speeds. Inter- 
estingly, a separate fit to the 3 decoupled spirals 
yields a bad match to all the end point locations 
expected from theory, that is, to the bar OLR, the 
arm CR, and the arm OLR. Since 2 of the objects 
(NGC 864 and NGC 3513), though, are consistent 
with ending at the arm CR, more data are needed 
to better understand the dynamics of spirals de- 
coupled from the bar. 

6. Conclusions 

Our results show that a connection exists be- 
tween bar/spiral dynamics and star formation. 
We have found indications of the existence of az- 
imuthal color (age) gradients across the bars and 
spirals of disk galaxies (although different mecha- 
nisms of star formation triggering may take place 
in both types of regions, see II. ip . Through the 
comparison of optical and NIR images with stellar 
population synthesis models, a link can be estab- 
lished between large-scale star formation in the 
disks and bar/spiral dynamics. 

For the bar regions, we compare the computed 
CR positions with the bar's end, and with results 
from other authors, both theoretical and obser- 
vational. The calculated CR radii for the bar 
pattern speeds are close to the bars' end points, 
in agreement with theoretical expectations. The 
analysis of azimuthal color (age) gradients shows 
that non-circular motions are important. In the 



^^This plot includes NGC 4314 and NGC 5383, although no 
gradients were found in the arms of these objects. 

•^^The bisector line was obtained by first fitting the 

OLS( Y|X) and OLS(X|Y) weighted by the er- 

rors |B evingto n fc RobinsonI l2003l ) . and then applying 
llsobe et al. I 1I1990I') formula for the OLS bisector slope. 



case of bar regions, the use of a circular dy- 
namic model produces a trend to overestimate 
ilp'^" for regions inside CR. This trend is similar 
to the one encountered in non-barred and weakly- 
barred spirals and, as already demonstrated by 
Martmez-Garcia et al. I (l2009bl) . does not imply the 
absence of a pattern speed. 

For regions in the spiral arms of barred galax- 
ies, we find that "inverse" color gradients (10 of 
20) also follow a trend that can be attributed to 
non-circular motions. In this case, though, the 
overestimation of £7^?"''' occurs near the CR ra- 
dius of the bar and converges to zero at higher 
radii. We also find gradients in the spiral arms 
where stellar aging follows the direction of rota- 
tion. The ripP"'*'' values derived from these regions 
are in general lower than fip'^''. 

Out of 9 galaxies with detected gradients in 
both the bar and the arms, 6 appear to have 
flp^" w fi^P"''*'; with one exception, these are all 
galaxies with early Hubble types. The remaining 
3 galaxies are late Hubble types, and appear to 
have rj);^' ^ n^P'"''^!. 

From the presence of azimuthal color (age) gra- 
dients alone, it is difficult to discern between mod- 
ern theories of the origin of spiral arms in barred 
galaxies; the pattern speeds that we can obtain 
based on the gradients, however, can provide sig- 
nificant clues. 
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Fig. 3. — NGC 718, region A. (a.) g-band deprojected mosaic of the galaxy, in logarithmic scale, (b.) Solid 
black line and left y-axis: observed Q{rJgi) profile vs. azimuthal distance, d, in kpc; dotted blue line and 
second- from- right y-axis: observed (g — J) color vs. d; dashed red line and rightmost y-axis: observed Kg 
sm'face brightness (mag arcsec"^) vs. d. (c.) Solid line: zoomed-in version of Q{rJgi) vs. d profile. Dotted 
line: stellar population model, "stretched" in tagc to fit the data; IMF Afuppcr = IQMq. Model stellar age, 
in units of 10^ yr, is shown in the upper a;-axis. Dashed line: model including "dissolution of stellar groups" , 
see ijS.l.ll Reduced values cover the same time interval for both models, (d.) Upper plot: isophote 
ellipticity vs. B}^"^; lower plot: isophotal PA vs. i?^/^. Hatched areas: bar corotation region (see text). 
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Fig. 4. — NGC 718, region B. (c.) Dotted and dashed lines: stellar population models, IMF Muppcr = IQMq. 
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Fig. 5. — NGC 718, region C. (c.) Dotted and dashed lines: stellar population models, IMF Muppcr = lOOAf©. 
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Fig. 6. — NGC 864, region A. (c.) Dotted and dashed lines: stellar population models, IMF Muppcr = lOM©. 
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Fig. 7. — NGC 864, region B. (c.) Dotted and dashed lines: stellar population models, IMF Muppcr = IOMq. 
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Fig. 9. — NGC 864, region D. (c.) Dotted and dashed lines: stellar population models, IMF Muppcr = lOM©. 
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Fig. 42. — NGC 3513, region D. (d): isophotes from _ff-band mosaic, (c.) Dotted and dashed lines: stellar 
population models, IMF Mupper = lOOM©. 
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Fig. 44. — NGC 3513, region F. (cf): isophotes from _ff-band mosaic, (c.) Dotted and dashed lines: stellar 
population models, IMF Mupper = lOM©. 
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is reached. During phase 2, the velocity dispersion increases with time and the "diffusion of stellar orbits" 
causes the dissipation of the stellar group. 



71 



1.7 - 



1.65 - 



L6 



L55 - 



1.5 - 



T ■ 1 r 



IMF upper mass limit: 100 Mg 
IMF upper mass limit: 10 M,^ 




J _L I I I L 



J I I I I I I I L 







5x10^ 10^ 

t(yr) 



1.5xl0« 



Fig. 47. — Same as figure l45l including dissolution of stellar groups (|Wielenlll977( ) after 50 Myr (see §[5TTT]) 



72 



1.8 - 



L7 - 



DO 



O 

I 



1.5 - 



1-4 - 



T — ^ — ^ — ^ 

/ 



T — r 



T — I — r 



T — r 



\ 




J I I L 







J L 



J L 



5x10^ 



t(yi-) 



1.5x108 



Fig. 48. — Index Q(r Jgi) vs. time for CB07 models, reddened as per the two-component dust model 
of , Chariot fc Falll (|2000t ). The duration of the burst is 2 x lO'^ yr, with a Salpeter IMF, 2% by mass of young 
stars, and solar metallicity. Lower and upper mass limits are 0.1 and 10 Mq, respectively. 



73 



36 - 




w 34 - 



20 



6 



8 10 
R (kpc) 



12 



14 



Fig. 49. — Pattern speed, ^Ip, for NGC 7479, obtained from the comparison between observations of color 
gradient candidates and stellar population synthesis models; the galactocentric radii of the color gradients 
are indicated in the x-axis. Regions A, B, and C belong to the bar, while regions D and E are located in the 
spiral arms. 
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Fig. 50. — Pattern speed, flp, for NGC 3513, obtained from the comparison between observations of color 
gradient candidates and stellar population synthesis models; the galactocentric radii of the color gradients 
are indicated in the x-axis. Regions A, B, and C belong to the bar, while regions D, E, and F are located in 
the spiral arms. 
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Fig. 53. — Same as figure [52l for the spiral regions only. Red open triangles: "inverse" color gradients; black 
solid circles: stellar aging and disk rotation have the same direction. 
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Fig. 54. — Mean arm maximum extent, i?^^™, obtained by eye from NIR data, vs. bar OLR, -Rqlr, 
assuming Rqk = l-2i?bar- Red open triangles: ri^^"^ ^ fip''™; black solid circles: ^^'^'^ « n^;™. Dotted line: 
^cnd — ^OLR' dashed line: OLS (weighted by errors) bisector line. 
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